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Molecular beam epitaxy of high structural quality Bi 2 Se3 on lattice matched 
InP(lll) substrates 
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Epitaxial layers of the topological insulator Bi2Se3 have been grown by molecular beam epitaxy on laterally 
lattice-matched InP(lll)B substrates. High resolution X-ray diffraction shows a significant improvement of 
Bi2Se3 crystal quality compared to layers deposited on other substrates. The measured full width at half 
maximum of the rocking curve is Aw = 13 arcsec, and the (u> — 29) scans exhibit clear layer thickness fringes. 
Atomic force microscope images show triangular twin domains with sizes increasing with layer thickness. The 
structural quality of the domains is confirmed on the microscopic level by transmission electron microscopy. 
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The V2VI3 compounds (such as Bi2Se 3 ) have been 
well known as thermoelectric materials for over 50 years. 
More recently, 3-dimensional (3D) topological insulator 
(TI) properties have been predicted for a number of 
these materials and the existence of a topologically pro- 
tected surface state has been confirmed by angle-resolved 
photoemission spectroscopy (ARPES)i— While the ex- 
istence of this surface state is well established, achieving 
a high structural quality crystal, free of any defects, is 
very challenging. In an attempt to reduce defect den- 
sity and to produce thin layers useful for devices, growth 
of Bi2Se3 layers by molecular beam epitaxy (MBE) has 
been extensively studied on various substrates, including 
Si(lll)^ GaAs(lll)^ Al 2 O 3 (110)£ SrTi0 3 (Hl)^ and 
CdS(0001)3 While the structural quality of these lay- 
ers has steadily improved, the narrowest rocking curve, 
an important indicator of structural quality, which has 
been reported to date, had a full width at half maximum 
(FWHM) of Alu ~ 360 arcsec^ i.e., below the quality usu- 
ally obtained in zincblendc semiconductor MBE layers.™ 
Here, we make use of nearly lattice matched InP(lll)B 
substrates with low miscut (< 0.2°). These have a hexag- 
onal surface lattice with a lattice-mismatch of only 0.16% 
to the ab-plane of Bi2Se3 i 12 ' 13 We show that these sub- 
strates allow for MBE growth of Bi2Se3 layers with a high 
crystal quality, which further increases with layer thick- 
ness. For a 250 nm thick layer, a high resolution rocking 
curve with FWHM of Auj ~ 13 arcsec is obtained. 

Bi2Se3 has a crystal symmetry belonging to the 
space group R3m (D| d ) and has lattice constants of 
a = b = 4.143 A and c = 28.636 AM It consists of so- 
called quintuple layers (QLs) Se-Bi-Se-Bi-Se with ABC 
stacking order along the [001] direction. The bonding 
type within a QL is covalcnt, in contrast to the weak van 
der Waals bonding between QLs.— 

Growth of the Bi2Se3 is carried out in a CreaTec UHV 
MBE system with a base pressure of 10 -10 mbar. Ele- 
mental Bi (6N) and Se (6N) are evaporated from standard 



effusion cells. The undoped InP(lll)B substrates have 
a miscut to the (111) plane specified to be below 0.2°. 
Atomic force microscope (AFM) measurements of a ref- 
erence substrate reveal terraces of about 1 /im length, 
corresponding to a miscut of 0.02°. Before growth, the 
substrate surface oxide is removed in 50% hydrofloric acid 
(HF). The substrate is then transported in an N2 atmo- 
sphere and quickly loaded into the MBE system. The 
growth is monitored by reflection high-energy electron 
diffraction (RHEED) using a CCD camera. After growth, 
the samples are characterized by a DME DualScope 95- 
50 atomic force microscope (AFM), a Philips X'Pcrt 
MRD diffractomcter for high resolution X-ray diffrac- 
tion (HRXRD) (the resolution in w-dircction is limited 
to about 6 arcsec by the instrumental broadening of the 
4xGe(220) monochromator and the acceptance angle of 
the analyzer crystal is 12 arcsec) and a FEI Titan 80-300 
scanning/transmission electron microscope (S/TEM) op- 
erated at 300kV. TEM cross-sectional specimens are pre- 
pared by focused ion beam milling. 

Previous experience has shown 300 °C to be the op- 
timal substrate temperature for growing Bi 2 Se3 in our 
MBE chamber. To reduce phosphorus out-diffusion from 
the InP substrate^ we choose a two-temperature growth 
start (2T-start)^£ to reach this temperature. We first 
heat the substrate to a temperature T s = 250 °C and 
cover the InP with approximately two QLs of Bi2Se 3 . 
The sample is then heated to T s = 300 °C and this ini- 
tial Bi2Se3 layer is annealed in Se atmosphere to im- 
prove its crystal quality. This is followed by deposi- 
tion of the bulk of the layer in normal MBE growth 
mode. As a control, we have also prepared a sample 
with a one-temperature growth start (lT-start), where 
we heat directly to T s = 300 °C and start growth. For 
both methods, the Bi2Se3 layer is grown at a rate of ap- 
proximately 1 nm per minute (0.17 As -1 ), using beam 
equivalent pressures (BEP) of ps e = 6 x 10 -6 mbar and 
PBi = 2 x 10~ 7 mbar. After growth, the samples are 
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cooled to T s = 140 °C in a Sc atmosphere to avoid for- 
mation of Se vacancies. Below, structural data on four 
samples, A, B, C, and D of different Bi2Se3 layer thick- 
nesses, corresponding to growth times (after anneal) of 
22 min, 1 h, 3 h, and 6 h are presented for the 2T-start 
procedure. Sample E is grown by lT-start procedure 
with a growth time of 1 h. 

RHEED is used to monitor the growth. Fig.[TJa) shows 
the intensity of the specular spot as a function of time. 
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FIG. 1. (a) RHEED specular spot intensity vs time of the 
Bi2Se3 growth start (T s = 250 °C) as well as the beginning 
of the second growth step (T s = 300 °C) for sample B, (b) 
and (c) RHEED images of the Bi2Se3 surface observed along 
the [011] and [121] directions of the substrate, after growth is 
completed, (d) - (g) AFM images (7x7 /im, lower right cor- 
ner: z scale bar labeled by the maximum height distribution 
in the image) of samples A, B, C, and D, respectively. 



We observe two oscillations assigned to the growth of 
about 2 QLs during the low temperature growth step. A 
slight increase of intensity continues when closing the Bi 
shutter, and saturates 25 s later. The sample is then an- 
nealed for - 1 h at T s = 300 °C, and normal MBE growth 
is initiated. At this high-temperature growth start dis- 
tinct oscillations are again observed. These gradually 
decrease in amplitude within the first 4 QLs. Extrap- 
olating the layer thickness from the period of these os- 
cillations and the growth time gives a film thickness of 
d_4 = 19 nm, d^ = 54 nm, Ac = 166 nm, djj = 315 nm, 
and d^ = 53 nm for samples A through E, respectively. 

RHEED streaks and Kikuchi finest are visible during 
the entire growth, and the FWHM of the streaks decrease 
with increasing layer thickness. The RHEED patterns in 
Figs.QJb) andQJc) are obtained along the [011] and [121] 
directions of InP, respectively, after the second growth 
step of the Bi2Se3 layer of sample B. The change in streak 
separation between the [011] and [121] direction is a fac- 
tor of v3, confirming the c-axis of the Bi2Se3 hexagonal 
unit cell is normal to the InP (111) planes. 

The AFM images of the samples (Figs.(Tfd)-(g)) reveal 
a surface comprised of domains of pyramidal shape, ro- 
tated by 60° with respect to each other, likely the result 
of twinning during heteroepitaxy. The valleys between 
these domains have a depth of a few tens of nm in the 
case of the 250 nm thick sample D. Distinct steps are 
observed on the flanks of the domains. Their height is 
~ 1 nm, which corresponds to single QL steps. The av- 
erage size of the triangles increases with layer thickness. 
This is presumably caused by domains of one orientation 
overgrowing domains of the other type, and thus, merg- 
ing into one larger domain. 

The (ui — 26) HRXRD scan of sample D presented 
in the inset of Fig. HJa) shows sharp and intense (00Z) 
(I = 3n, n € Ni) peaks of Bi2Sc3 in addition to the 
(111) and (222) peaks of the substrate. The absence of 
other peaks confirms the growth of primarily single stoi- 
chiometric phase Bi2Se3. The higher resolution (cj — 29) 
scans of the (003) reflection (shown in Fig. &1&)) reveal 
pronounced layer thickness fringes. Their period corre- 
sponds to a layer thicknesses of Aa = 24 nm, d# = 63 nm, 
dc = 159 nm, d^> = 250 nm, and d^ = 66 nm for samples 
A, B, C, D, and E, respectively. These are not inconsis- 
tent with the thicknesses estimated from RHEED when 
taking into account that the latter assumes a perfectly 
constant growth rate. 

A HRXRD rocking curve (w-scan) of sample D 
shows a very sharp (003) Bi2Se3 peak with a FWHM 
Auj d = 13 arcscc (Fig.gfb)). This FWHM is compara- 
ble to the instrumental resolution. The inset of Fig. [2jb) 
shows the measured FWHM Au>, which decreases with 
increasing layer thickness. This shows that the mosaic- 
ity of the layer decreases with thickness. The FWHM of 
samples E and B with similar layer thickness, but with 
different growth starts, is comparable. This suggests that 
not the annealing of the starting layer in the 2T-start, of- 
ten used for improving the crystal quality on lattice mis- 
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matched substrates^ but rather the InP(lll) substrate 
results in the achieved high structural quality. The sharp 
(00 I) peaks are surrounded by a weak diffusely scattered 
background and shoulders in w-direction. The data are 
fitted by two pseudo-Voigt profiles with 14% integrated 
intensity for the background. 

The reciprocal space map in Fig. Eta) clearly shows the 
narrow (003) reflection of the 250 nm Bi2Se3 layer (sam- 
ple D) along with the associated layer-thickness fringes 
and the shoulders, as well as a diffuse background feature, 
which is elongated in the Qh direction. The background 
and shoulders could be related to any combination of fi- 
nite size effects^ of the domains, tilts at domain bound- 
aries, interface defects, or the slight "waviness" of the QLs 
seen in TEM images (see Fig. 0Jb)). We note, however, 
that a direct comparison of our HRXRD data with the 
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FIG. 2. HRXRD of the Bi 2 Se 3 (003) reflection on a loga- 
rithmic scale, (a) (w — 29) diffractograms for layers of var- 
ious thicknesses, showing distinct thickness fringes (subse- 
quent curves offset by two decades for clarity). The inset 
shows a large-scale (cj — 26) scan with the peaks identified, 
(b) rocking curve scan of the 250 nm sample D with a FWHM 
of 13 arcsec, the profile fit (blue) is the superposition of two 
pseudo-Voigt fits (red) , the inset shows the FWHM of rocking 
curves of the various samples plotted versus layer thickness. 
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FIG. 3. (a) Reciprocal space map of the Bi2Se3 (003) reflec- 
tion of sample D. (b) Asymmetric reflections of {0 1 5} layer 
planes and {0 2} substrate planes (curve offset by three 
decades for clarity) as a function of in-plane rotation cp, the 
inset shows a close up of a layer peak. 



published data on epitaxial Bi2Sc3 layers is difficult, as 
the reported w-widths may be increased by instrumental 
resolution. 

Asymmetric reflections from the {0 1 5} Bi2Se3 layer 
planes and the {0 2} substrate planes are measured 
as a function of in-plane rotation <f>, in order to study 
twinning and the influence of the substrate's lattice on 
the layer. The Bi 2 Se3 reflections occur every 60° (see Fig. 
|3{b)), instead of the 120° expected for trigonal symmetry, 
consistent with reports of hot wall epitaxy of Bi2Se3 on 
InP(lll)BJ^ The observed six fold symmetry is caused 
by twin domains with either stacking order ABC or the 
reversed ACB. The three {0 2} reflections of the sub- 
strate are at the same azimuthal angle cf> as those of one 
triplet of the Bi2Se3 layer. This confirms that the hexag- 
onal lattice of the Bi2Se3 layer is oriented parallel to that 
of the substrate, to within a twist of order 0.07°, given 
by the FWHM of the Bi2Se3 reflections in 4> direction 
(see inset Fig. 03h)). Apparently, the van der Waals 
bonds between the Se and the substrate are sufficient to 
orient the hexagonal lattice structure, but do not deter- 



4 





FIG. 4. (a) Dark-field TEM image of sample B. White arrows 
indicate the positions of twin domain boundaries, (b) cross- 
sectional high-angle annular dark field S/TEM image of the 
Bi2Se3/InP interface region within a single domain. 



mine the choice of "B" or "C" sites for the second atomic 
layer (Bi), allowing both ABC and ACB stacking, and 
thus the formation of twin domains. 

TEM studies of sample B further characterize layer 
and interface properties. Selected area electron diffrac- 
tion patterns display sharp spots and confirm the parallel 
orientation of the Bi2Se3(001) and InP(lll) planes. Both 
electron diffraction patterns and dark-field TEM images 
reveal the presence of rotational and lamellar twin do- 
mains in the layer (Fig.[4|a)).— The observed decreasing 
number of twin domain boundaries (indicated by white 
arrows) with distance to the interface supports the pre- 
viously suggested mechanism of merging and overgrowth 
of twin domains, which results in a domain size increas- 
ing with layer thickness. At the bottom interface, there 
is a layer of relatively low crystalline quality along with 
higher crystalline quality regions up to a QL thickness, 
which is followed by well-crystallized QLs (Fig. BJb)). 

In summary, using InP(lll) substrates, which arc 
nearly lattice-matched to the natural lattice constant 
of Bi2Sc3 allows for the MBE growth of layers of high 
structural quality. HRXRD investigations reveal that 
the layers have low mosaicity-tilt, and -twist, and are 
of uniform thickness across the sample. AFM and 
TEM studies show that the layers are nevertheless not 
truly monocrystalline but rather comprised of pyramid 
shaped twin domains of a size which increases with layer 
thickness. The formation of twin domains appears to 
be an intrinsic property of heteroepitaxy of Bi 2 Se 3 on 
hexagonal surfaces like InP(lll) due to the two possible 
orders of layer stacking in the Bi 2 Se3 unit cell, i.e., 
ABC ABC... and ACB ACB.. H It would be of interest 



to study growth on InP substrates of various orientations 
to see if this property can be controlled. Finally, the 
length scale of the achieved domains is already sufficient 
to consider patterning mesoscopic devices such that they 
comprise of a single high quality domain. 
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